Structural Transitions and Global Minima of Sodium Chloride Clusters 
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In recent experiments on sodium chloride clusters structural transitions between nanocrystals 
with different cuboidal shapes were detected. Here we determine reaction pathways between the low 
energy isomers of one of these clusters, (NaCl)35Cl~. The key process in these structural transitions 
is a highly cooperative rearrangement in which two parts of the nanocrystal slip past one another on 
a {110} plane in a (110) direction. In this way the nanocrystals can plastically deform, in contrast 
to the brittle behaviour of bulk sodium chloride crystals at the same temperatures; the nanocrystals 
have mechanical properties which are a unique feature of their finite size. We also report and 
compare the global potential energy minima for (NaCl^Cl - using two empirical potentials, and 
comment on the effect of polarization. 
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I. INTRODUCTION 



There has been much interest in phase-like transitions 
in clusters. Most of this research has been directed at the 
finite analogue of the solid- liquid transition [jp , whereas 
the possibility of transitions between different solid forms 
has received less attention [^pQ. Some examples are 
known: when supported metal clusters are observed by 
electron microscopy, their structure can change between 
icosahedral, decahedral and close-packed ||||. However, 
the role of the electron beam and the surface is uncer- 
tain. Similar transitions have also been suggested for 
free atomic clusters of certain sizes, which are driven by 
entropy differences between the icosahedral, decahedral 
and close-packed morphologies 0-^). However, it is not 
clear how these changes could be detected even if the 
large free energy barriers involved |l0| are surmountable. 

Coexistence between solid-like isomers of small bi- 
nary salt clusters has previously been described by Rose 
and Berry More recently, clear examples of struc- 
tural transitions have emerged from experiments on NaCl 
clusters. These clusters have only one energetically 
favourable morphology: the magic numbers that ap- 
pear in mass spectra correspond to cuboidal fragments 
of the bulk crystal (rocksalt) lattice [pd|-p^[, hence the 
term nanocrystals. Indirect structural information comes 
from the experiments of Jarrold and coworkers which 
probe the mobility of size-selected cluster ions. For most 
(NaCl)ArCl~ with N > 30, multiple isomers were de- 
tected which were assigned as nanocrystals with different 
cuboidal shapes [H . The populations in the different iso- 
mers were not initially equilibrated, but slowly evolved, 
allowing rates and activation energies for the structural 
transitions between the nanocrystals to be obtained [15| . 
Based on the small values of the activation energies, Hud- 



gins et al. suggested that the rearrangement mechanisms 
might involve a sequence of surface diffusion steps. 

The aim of the present work is to examine this hypoth- 
esis by identifying the mechanisms for one of the clusters 
that was studied experimentally, (NaCl^Cl - . We also 
report the global potential energy minima of (NaCl) atC1~ 
for two different empirical potentials which enables us to 
comment on the effects of including polarizabilities and 
hence induction energies. Methods are outlined in section 
U and the global minima are discussed in HI . Then in sec- 
tion ^ we describe the mechanism that we found to me- 
diate the structural transitions in (NaCl) 35 CI - . Finally, 
in section ^ we discuss the implications of this mecha- 
nism for the mechanical properties of NaCl nanocrystals 
and its relevance to other alkali halides. 



II. METHODS 
A. Searching the PES 

The identification of mechanisms for the structural 
transitions of (NaCl) 35 CP presents a considerable chal- 
lenge to the theoretician since the half-lives of the least 
stable isomers are of the order of milliseconds |jl5| , 
whereas the time scales that can be probed by conven- 
tional molecular dynamics simulation are only on the 
order of nanoseconds. The difficulty is that in molec- 
ular dynamics or Monte Carlo simulations the system 
spends most of the time vibrating about a minimum on 
the PES with transitions between minima occurring only 
rarely. One approach to enhancing the rate of occurrence 
of rare events, such as these structural transitions, is to 
bias the system towards the transition regions using um- 
brella sampling p6[; this method is particularly suited 
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to calculating free energy barriers. Here we use an ap- 
proach which is ideal for finding reaction pathways for 
complex processes. Previously it has been used to find 
a pathway between the face-centred cubic global mini- 
mum and the lowest energy icosahedral minimum of a 
38-atom Lennard- Jones cluster and to identify re- 
laxation mechanisms in amorphous silicon E^ j. Steps 
are taken directly between minima [ |i7||l^ , |2C| 7 thus al- 
lowing large distances on the PES to be traversed. To 
do this, we first find a transition state connected to the 
current minimum using the eigenvector-following tech- 
nique pl|-|2^| . We then calculate the corresponding rear- 
rangement mechanism and thereby obtain the new min- 
imum. Finally, we decide whether to accept the step to 
the new minimum, typically on the basis of a Metropo- 
lis criterion p7f| . By repeating this process the system 
performs a walk amongst connected minima on the PES. 

The steps between adjacent minima are not biased to 
go in any particular direction. Nevertheless, by perform- 
ing an extensive search of the low-lying minima we were 
able to find paths between the various rocksalt-type iso- 
mers of (NaCl) 35 Cl~. In the process over 4500 minima 
and 5500 transition states were characterized. Although 
there is no guarantee that we have found the shortest or 
lowest barrier rearrangements between the nanocrystals, 
we are confident that our paths are good estimates and 
representative of the optimal paths. 
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(There is a small typographical er- 
ror in equation (3) of reference The above formu- 
lae are given in atomic units. Adopting the notation of 
Stone [ pT| , component a of the induced dipole vector at 
site B due to site A, fj, B , is: 



a B ,F A ,(B) 



where A and B are the position vectors of the respective 
sites. The total induced moment in the present case is: 
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and 



R = IRI 



Two approaches are available for finding the dipole mo- 
ments, namely iterating the equations for fi B to self- 
consistency or rearranging the equation as 

M/x = Y so that (j, = M _1 Y, 



B. Potentials and geometry optimization 

Two popular empirical potentials for NaCl were con- 
sidered in the present work. The first is the Tosi- 
Fumi parameterization of the Coulomb plus Born-Meyer 
(C+BM) potential & 



Aije 



-ni/p 



where qi is the charge on ion i, is the distance be- 
tween ions i and j and Aij and p and parameters |p8| . 
We have also considered the more complex potential fit- 
ted by Welch et al. |29] for which the full vector form has 



previously been given by Phillips et al. 1 30 1 : 
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where 



Qi Qj 



where fi = ([Jr,fx , . . .). In the present work matrix in- 
version was used to obtain the self-consistent dipole mo- 
ments and the first analytic derivatives of the energy. 

The inclusion of polarizabilities in the Welch poten- 
tial makes this functional form much more expensive to 
evaluate than the simple C+BM form. Hence we con- 
ducted the most extensive searches of the PES with the 
latter potential and then reoptimized stationary points 
and pathways with the Welch form. Global minima were 
located using a guiding function approach for the Welch 
potential, as discussed in the next section. Transition 
states were located using a modified eigenvector-following 
approach. The basic algorithm has been described be- 
fore |^,^|, and was used in the present work with nu- 
merical second derivatives for the Welch potential. We 
also employed a new approach which does not require 
second derivatives and is more efficient [B4f. 



III. GLOBAL OPTIMIZATION 



Global potential energy minima were located for 
(NaCl^Cl - up to N — 35 for both of the empirical 
potentials described above. We employed the 'basin- 
hopping' or Monte Carlo minimization f35} | technique 
which has recently been investigated for a variety of 
atomic and molecular clusters HBQ-fill. 
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FIG. 2. Global minima for the Welch potential corresponding to local rather than global minima of the Coulomb plus 
Born-Meyer potential. Lines have been drawn using a distance cut-off to guide the eye. The sodium ions are represented by 
the smaller, darker circles and the chloride ions by the larger, lighter circles. 



Canonical Monte Carlo (MC) sampling was used to ex- 
plore the transformed PES, as described elsewhere P,p6|. 
For the C+BM potential five runs of 5000 MC steps 
were performed for every cluster size with each run start- 
ing from a random point. Two short runs of 200 steps 
each were also initiated using the lowest minima for 
(NaCl) 7V _ 1 Cl~ and (NaCl) Af+1 Cl~ as seeds. The maxi- 
mum step size for the displacement of any Cartesian co- 
ordinate was dynamically adjusted to give an acceptance 
ratio of 0.5 for a temperature corresponding to 316 K 
(474 K for the larger clusters). Final values for the max- 
imum displacement were typically around 1.5 A. To re- 
strict the configuration space to bound clusters we reset 
the coordinates to those of the current minimum in the 
Markov chain at each step . 

Since the Coulomb potential is long-ranged these clus- 
ters actually represent a rather easy global optimization 
problem because there are fewer minima on their poten- 
tial energy surfaces 0,^,^3) . Hence we are confident that 
most of the lowest energy structures we have identified 
for the simpler potential are the true global minima. In 
fact, the lowest minimum was usually the same for each 
of the five runs at every size considered. The Welch po- 
tential is much more time consuming to evaluate, and 
in this case only three runs of 4000 steps each were per- 
formed. We also employed a guiding function technique, 
as suggested by Hartke E3], where the simpler C+BM 
was used for partial geometry optimization followed by 
relaxation with the full Welch potential in every quench. 
Once again the same lowest energy structures were usu- 



ally found in each of the three runs for every size. The 
results also agree with calculations performed without a 
guiding function up to N — 17, and with those of Phillips 
et al. HH up to N — 14, although our energies are sys- 
tematically lower than theirs, presumably due to differ- 
ent unit conversion factors. Our energies and geometries 
are all well converged with the root-mean-square force re- 
duced to less than 10~ 8 hartree/bohr for every minimum. 
The most difficult case in this size range appears to be 
N = 31 where the 7x3x3 rocksalt global minimum 
was only found in two out of three runs for the Welch 
potential. 

All the results will be provided in a downloadable for- 
mat from the Cambridge Cluster Database fig] . The 
energies of the lowest minima found for both the C+BM 
and Welch potentials are given in Table | and the struc- 
tures are illustrated in Figures [l] and ^. The lowest ten or 
so local minima found for the B+CM potential were also 
relaxed separately under the Welch potential to establish 
the correspondence between local minima. 

For most sizes in the present range the structure of 
the global potential minimum is the same for the two 
potentials. However, the global minimum of one poten- 
tial is only a local minimum for the other potential at 
N = 6, 10, 11, 15, 18 - 20, 23, 24 and 29. For N = 2 
the global minimum is linear for C+BM but somewhat 
bent for Welch. For N — 5 the Welch global minimum 
has lower symmetry than the C+BM structure; the ap- 
parent D^d symmetry is actually slightly broken on close 
inspection. The similarity of the results for the two po- 
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tentials shows that the C+BM form, despite its simplic- 
ity, can give a reasonably reliable guide to the structure 
of these clusters and provides justification for our use of 
the C+BM potential to perform surveys of the energy 
landscape. 

As expected, the lowest energy structures have pre- 
dominantly rock-salt structures, and the particularly sta- 
ble sizes occur when complete cuboids can be formed 
(Figure |). These sizes (N = 4,13,22,31) agree with 
the magic numbers observed in the mass spectral abun- 
dance distributions p3| . It is interesting to note that 
at a number of sizes when complete cuboids cannot 
be formed a column of hexagonal rings appears (e.g. 
N = 16, 20, 25, 27). One of the more unusual structures 
is the (NaCl)24Cl~ global minimum for the Welch poten- 
tial. It has threefold symmetry with a trigonal bipyramid 
in the middle containing three (equatorial) Na + ions and 
two (axial) CI - ions. 




FIG. 3. The energies of the global minima as a func- 
tion of size for the Welch (solid line) and C+BM (dashed 
line) potentials. To make the size dependence more clear, 
we have used a function, E nve = a + bN 1/3 + cN 2/s + dN, as 
the energy zero. The coefficients of -Eavo have been chosen 
to give the best fit to the energies. For the Welch potential 
a = -0.5014, b = 0.2139, c = 0.0926 and d = -7.6829. For 
the C+BM potential a = -1.6782, b = 1.8878, c = -0.2828 
and d = -7.6360. 



IV. REARRANGEMENTS OF (NaCl) 35 Cr 

In the experiments on (NaCl^sCI - three peaks were 
resolved in the arrival time distribution. These peaks 
were assigned on the basis of calculated mobilities as an 
incomplete 5x5x3 cuboid, an incomplete 5x4x4 cuboid, 
and an 8 x 3 x 3 cuboid with a single vacancy Jl4] , p^[ . 
However, the lowest energy minima that we found for 
this cluster could be divided into four types (Table [l| 
and Figure ^) , namely the above three nanocrystals and 
a 6 x 4 x 3 cuboid with a single vacancy. As the latter 



structure was not considered when the structural assign- 
ments were made, its mobility has now been calculated 
using the exact hard-spheres scattering model |4(||^J . 

With these new data better agreement between the 
calculated and observed mobilities can be obtained by 
assigning the three experimental peaks to the 6x4x3, 
5x5x3 and 8x3x3 nanocrystals in order of increasing 
drift time j47|. Furthermore, this reassignment gives bet- 
ter agreement with the energetics. In the experiment the 
clusters convert to (what we now assign as) the 5x5x3 
nanocrystal as time progresses, indicating that this struc- 
ture has the lowest free energy. In our calculations the 
global potential energy minimum is also a 5 x 5 x 3 iso- 
mer. Moreover, the 6x4x3 nanocrystal is second lowest 
in energy; it is 0.06 eV (0.10 eV for the Welch potential) 
above the global minimum, whereas the corresponding 
value for the lowest energy 5x4x4 isomer is 0.21 eV 
(0.25 eV). 

The disconnectivity graph (or tree) Q,fl9| in Figure || 
provides a graphical representation of the PES. It shows 
which minima are connected by pathways below any 
given total energy. The end of each line represents a 
minimum on the PES and each node occurs at the energy 
of the transition state which first connects the two (sets 
of) minima. The tree separates the minima according to 
the particular rocksalt cuboid quite cleanly. This result 
shows that the barriers between minima with the same 
basic cuboidal shape are generally smaller than those be- 
tween the different types of nanocrystal. This separation 
holds least well for the 5x4x4 nanocrystals because of 
the many different ways that the nine vacant sites can be 
arranged. For example, the two lowest energy 5x4x4 
minima, L and O, have very different structures (Figure 
|]) and the energy at which the set of minima associated 
with O becomes connected to the 6x4x3 minima is 
lower than the energy at which it becomes connected to 
the set of minima associated with L. Minimum O is also 
much closer in configuration space to the lowest energy 
6x4x3 minimum than it is to minimum L (Table |l|) . 
Another example is provided by the two 5x5x3 iso- 
mers on the left of the figure which are separated by a 
very large barrier from the rest of the 5x5x3 structures 
because the ions occupy the opposite lattice sites to the 
other minima (i.e. the sodium ions, not the chloride ions, 
are located at the corners of the nanocrystal). 

One helpful way to characterize the topography of an 
energy landscape that has come from the protein folding 
community is in terms of funnels |^,^l|. A funnel is a 
set of pathways that converge to a low energy minimum. 
It has been suggested that a single deep funnel underlies 
the ability of proteins to fold to a unique native structure. 
Figure || shows that on the (NaCl^Cl - PES there are 
separate funnels corresponding to the 5x5x3, 6x4x3 
and 8x3x3 nanocrystals and a number of small funnels 
associated with the 5x4x4 nanocrystal. 
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FIG. 4. The lowest energy (NaCl) 35 Cl~ minima of the four types of nanocrystal for the Coulomb plus Born-Meyer (A, D, I 
and L) and Welch (A, E, I and O) potentials. D and E differ only in the position of the vacancy at an Na + site. The sodium 
ions are represented by the smaller, darker circles and the chloride ions by the larger, lighter circles. The minima are labelled 
as in Table [fi|. 




FIG. 5. Disconnectivity graph showing the multiple-funnel character of the (NaCl)3sCl~ PES. The branches that end at the 
100 lowest energy minima are shown. They are marked by the type of structure: 5x5x3: diamonds, 6x4x3: squares, 8x3x3: 
crosses, and 5x4x4: unlabelled. The dashed lines approximately divide the tree into the different types of nanocrystal. The 
numbers denote the number of minima in a branch. We also label the six minima that appear in Figure ^ with the appropriate 
letter. The energies are for the Coulomb plus Born-Meyer potential. 
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FIG. 6. Energy profiles of pathways between different 
(NaCl) 35 Cl~ isomers, (a) The shortest path and (b) the 
lowest energy path between the lowest energy 5x5x3 and 
6x4x3 isomers, (c) Two pathways between the lowest 
energy 6x4x3 and 8x3x3 isomers; the shorter pathway 
is the pathway with the lowest barrier. The sections of the 
profile with solid lines are illustrated in Figures fj] and ^| The 
energies are for the Coulomb plus Born-Meyer potential. 



This characterization of the PES helps to explain the 
initial presence of the metastable isomers in the exper- 



iments by Jarrold and coworkers. [ pL 4|jl5[ | The energy 
landscape efficiently 'funnels' high energy structures into 
rocksalt-type minima |32]]5|]. However, the barriers be- 
tween the different funnels are large, leading to a separa- 
tion of time scales for relaxation down the PES to a rock- 
salt structure and conversion of the metastable nanocrys- 
tals to the global minimum. In Ref. |54| we predicted that 
this behaviour would occur for alkali halide clusters. 

To understand the mechanisms for the structural tran- 
sitions we have visualized many pathways connecting 
nanocrystals with different dimensions. In virtually all 
the pathways the major shape changes are achieved by 
the same type of mechanism. A typical example of this 
process occurs in the shortest path between the lowest 
energy 5x5x3 and 6x4x3 isomers. The path involves 
four sequential transition states (Fig. |]a), the first three 
of which mediate local rearrangements of the empty sites 
in the 5x5x3 cuboid. The fourth rearrangement causes 
the major shape change and is depicted in Fig. 0a. In this 
'glide' mechanism the two halves of the cluster slide past 
one another on a {110} plane in a (110) direction. Other 
examples of this mechanism are illustrated in Figures 03 
and 0. 

Although the path shown in Figure 0a and 0a is the 
shortest between the lowest energy 5x5x3 and 6x4x3 
minima it is not the lowest in energy. The latter pathway 
(Figure 0b and 0) is more complicated; rather than pass- 
ing directly between the two nanocrystals it goes via a 
number of 5 x 4 x 4 minima, some of which are structurally 
similar to minimum O. Again the glide mechanism me- 
diates the main structural changes between a 5 x 5 x 3 
and a 5 x 4 x 4 nanocrystal (top of Figure 0) and from a 
5x4x4 to a6x4x3 nanocrystal (bottom of Figure 0). 
The former mechanism is not quite so clear-cut since the 
glide is coupled with local defect motion on the far side 
of the crystal. The middle rearrangement in Figure ac- 
tually involves the highest energy transition state on the 
pathway. In this rearrangement a triangle of ions slides 
over the surface of the nanocrystal. The final rearrange- 
ment in the pathway (not depicted in Figure 0) is the 
conversion of minimum E to D by the motion of a single 
sodium ion. 

The lowest energy pathways to the 8x3x3 nanocrys- 
tal all occur via the same transition state The pathway 
from the 6x4x3 nanocrystal is illustrated in Figures 
0c and 0b. Two glide rearrangements convert the lowest 
energy 6x4x3 isomer into an 8 x 3 x 3 isomer. In the 
final rearrangement on this pathway a single sodium ion 
moves so that the vacancy in the 8x3x3 nanocrystal 
occupies the site of lowest energy. 

It is not hard to see why the glide mechanism is 
favourable. On either side of the {110} plane are rows 
of oppositely-charged ions (these can be clearly seen on 
the top face of the nanocrystals in Fig. 0a and b) . When 
the nanocrystal slides on this plane in the direction of 
these rows, no ion comes any closer to the nearest ion 
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of the same charge. This situation only holds for {110} 
planes and the (llO) direction. The activation barrier 
arises from the loss of some favourable contacts between 
oppositely-charged ions. It is also interesting to note that 
this plane and direction correspond to the primary slip 
system for dislocation motion in bulk NaCl crystals [p5j . 

One other mechanism by which major shape changes 
can occur is illustrated in Figure fflc. In this 'hinge' pro- 



cess the two halves of the cluster split apart and rotate 
about a common edge until they meet up again. Al- 
though spectacular this mechanism is unlikely to be rel- 
evant to the experiments because the barrier is probably 
too large (Fig. ||b) to be surmountable at the appropriate 
temperature; the transition state is 1.57eV (1.21 eV for 
the Welch potential) above the 8x3x3 minimum. 



(a) 



5x5x3 



6x4x3 








FIG. 7. Portions of the three rearrangement mechanisms corresponding to Fig. |fja and c. (a) A glide rearrangement which 
converts a 5 x 5 x 3 into a 6 x 4 x 3 isomer, (b) Two glide rearrangements which are part of the lowest energy path between 
the 6x4x3 and 8x3x3 isomers, (c) A higher energy 'hinge' mechanism which converts a 6 x 4 x 3 isomer directly to an 
8x3x3 isomer. 



From the experiments of Jarrold and coworkers rates 
and activation energies were obtained for the conver- 
sion of the metastable nanocrystals to the global min- 
imum fl5f| . Therefore, it would be useful if we could 
make estimates of the activation energies to provide in- 



dependent support for our results. However, this is non- 
trivial. There are a number of possible approaches. For 
a complex process one needs to consider not the energy 
barrier, but the free energy barrier. The latter barrier 
could be computed by umbrella sampling |l(| if a suit- 
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able order parameter was found which could distinguish 
the different nanocrystals. The rate of barrier crossing 
could then be calculated by performing a series of simu- 
lations which start from the top of the barrier |56||37j . A 
second approach is to use a master equation to model the 
flow of probability between the minima on the PES [p8j . 
However, to achieve reasonable results one would need 
a larger set of minima and transition states than we 
have obtained here. Furthermore, this approach also re- 
quires an expression for the rate constants for transitions 
between minima. Although these are easily calculated 
within the harmonic approximation using RRKM theory, 
more accurate rate constants which include the effects of 
anharmonicity are far harder to obtain. There are also 
two recently-developed techniques which could possibly 
be applied to calculate the rate constants: the macrostate 
variational method |p9| and the transition path sampling 



approach of Chandler and coworkers |6y,|Gl| . 

All these approaches are computationally demanding 
and beyond the scope of the present work. However, we 
can obtain an estimate of the activation energy if we make 
the approximation that the only important path between 
the nanocrystals is the one with lowest energy. The acti- 
vation energy can then be equated with the energy differ- 
ence between the highest energy transition state on this 
path and the starting minimum. At absolute zero this en- 
ergy barrier is equal to the free energy barrier, however 
at non-zero temperatures the free energy barrier is likely 
to be reduced by the entropy associated with the multi- 
plicity of paths between the two states. This interplay of 
energy and entropy has been observed in the free energy 
barrier between the lowest energy face-centred-cubic and 
icosahedral minima for a 38-atom Lennard- Jones clus- 
ter @. 



5x5x3 



5x4x4 






5x4x4 





6x4x3 




FIG. 8. Three rearrangements from the lowest energy pathway between the lowest energy 5x5x3 and 6x4x3 minima, 
for which the energy profile is shown in Fig. |^b. The middle frame is the transition state and on either side are the two minima. 
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In Table III we give our estimates for the activation en- 
ergies of pathways between the lowest energy isomers of 
each cuboidal nanocrystal. In the experiment it was only 
possible to find the activation energies for rearrangements 
between the 6x4x3 and 8x3x3 nanocrystals and the 
5x5x3 nanocrystal; the values were 0.57 ± 0.05 eV and 
0.53 ± 0.05 eV, respectively |l5). These values are of the 
same order of magnitude as our estimates of the activa- 



tion energies (Table [II), although somewhat lower than 
we find for these specific transformations. The differ- 
ences could be due to the empirical potentials employed 
or the approximations involved in our estimation of the 
activation energy, which would probably produce over- 
estimates, as observed. It is also possible that we may 
not have found the lowest energy pathways. Further- 
more, the barriers are calculated with respect to the low- 
est energy minima for a particular type of nanocrystal. 
However, as the barriers between minima with the same 
cuboidal shape can be quite large (e.g. in Figure ||a there 
is a barrier between two 5x5x3 isomers which is 0.66eV) 
it is not certain that the nanocrystals in the experiment 
are in the lowest energy minimum for that shape. 

Given these uncertainties it is not yet possible to use 
comparison of the experimental and estimated activation 
energies to provide additional confirmation of the glide 
mechanism. However, the universality of the glide mech- 
anism amongst the many low energy paths that we have 
discovered strongly suggests that this process does medi- 
ate the structural transitions. 

Although most of the results in this section are for the 
C+BM potential those we have obtained using the Welch 
potential encourage us to believe that the C+BM form 
provides a reasonable description of the interactions. The 
global minimum for (NaCl) 35 Cl~ is the same for the two 
potentials. Furthermore, the low energy minima are sim- 
ilar (Table ||. There are a few changes in the energetic 
ordering of the minima; most significantly the 5x4x4 
minimum O is significantly lower than minimum L for the 
Welch potential. It is also interesting to note that the es- 
timates of the activation energy are systematically lower 
when the Welch potential is used (Table |lll|) bringing the 
values closer to those found in experiment |15|. 



V. DISCUSSION 

At the transition state of the glide mechanism there 
are stacks of hexagonal rings perpendicular to the direc- 
tion of sliding (Figure ^) . Hexagonal rings are a common 
motif in small alkali halide clusters |52], and when the 
ratio of the cation to anion size is smaller than for NaCl, 
e.g. in the lithium halides, NaBr and Nal, the hexagonal 
structures are more stable than for NaCl |63| . Hence the 
glide mechanism might be even more favourable in these 
systems. To prove this supposition calculations similar 
to those in this paper could be performed. Furthermore, 



if this suggestion is correct the dependence of the kinet- 
ics of structural transformation on the particular alkali 
halide could be then used as an experimental test of the 
mechanism we have found here. 

As the glide mechanism is a cooperative mechanism, 
unlike the surface diffusion mechanism originally sug- 
gested [^5|, it is expected that the barriers for this type 
of rearrangement would generally increase with size (al- 
though as with any cluster property it is likely that 
specific size effects would also be superimposed on this 
trend). The structural transitions would become more 
difficult with increasing size. Again this could probably 
be used as an experimental test of the mechanism. 

The glide mechanism also has implications for the me- 
chanical properties of these nanocrystals. It allows them 
to deform spontaneously and plastically to a thermo- 
dynamically more stable structure; in other words, the 
nanocrystals are soft. Similar homogeneous slip mecha- 
nisms, where whole planes of atoms slide past one an- 
other, have been seen in simulations of strained metal 
nanowires (6jj6^]. As the barrier to these cooperative 
processes increases with the area of the sliding surface, 
they become less feasible as the size of the system in- 
creases. For the metal nanowires the increased barrier 
leads to a change in the slip mechanism to a more lo- 
calized process, namely dislocation motion, with increas- 
ing size [p5[; the dislocations maintain the ductility of 
the metal nanowire. In contrast, for NaCl nanocrystals, 
at temperatures relevant to the experiments considered 
here (7-67°C |l4j]), dislocation motion is much more dif- 
ficult than for metals. Therefore, the increasing barrier 
appears to lead to a dramatic change in mechanical prop- 
erties. The nanocrystals become harder as their size in- 
creases, until a point is reached where an applied stress 
is more likely to cause fracture than plastic deformation, 
and the typical brittle behaviour of bulk NaCl crystals is 
recovered. 

The mechanical properties of these NaCl nanocrystals 
provide another example of the unique finite-size prop- 
erties of clusters. It might be possible to confirm the 
increased plasticity we predict using a microscope tip to 
deform NaCl nanocrystals soft-landed on a surface. 

Previously, Ball et al. have compared relaxation to the 
global minimum of (KC1)32 with Ari9 in terms of 'mono- 
tonic sequences' |p2[ . On this basis the salt cluster was 
described as a 'structure-seeker' because rapid removal of 
kinetic energy usually leads to rocksalt structures; in con- 
trast it is relatively easy to quench Arig into a defective 
double icosahedron. The present results for (NaCl) 35 Cl~ 
are relevant to this relation between the energy landscape 
and relaxation dynamics, regardless of whether the glide 
mechanism is operative in (KC1)32- 

Assuming that the model (KC1)32 cluster discussed 
elsewhere has a landscape similar to (NaCl) 35 CP wc 
would expect different cuboidal rocksalt morphologies 
to lie at the bottom of separate funnels. Relaxation 
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within each funnel is efficient, but relaxation between 
funnels will probably occur on a much longer time scale. 
Hence we can call the salt clusters 'structure-seekers' if by 
'structure' we mean any of the cuboidal minima. By the 
same token, then, we should perhaps have bunched the 
defective double icosahedra together with the global min- 
imum for Arig in the previous study p2[ . In fact Arig is 
a 'magic number' cluster with an essentially single funnel 
landscape and so Arig is also quite an efficient 'structure- 
seeker'. A better contrast would be provided by a cluster 
bound by a short range potential, where the landscape 
would be rougher, or by Ar38, whose double funnel land- 
scape we have investigated elsewhere |^- |lo|Jl7| , ^9|] . The 
point is that the salt cluster relaxes faster down any one 
funnel than Arig, but may not reach the global minimum 
except over a longer time scale. 
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TABLE I. Energies of the lowest minima found for 
(NaC^jyCP clusters with the Coulomb plus Born- Meyer and 
Welch potentials. 



N 






welcn 




Energy /eV 


PG 


Energy /eV 


PG 


1 


-7.712 


Dooh 


-7.960 


Dooh 


2 


— 14.948 


Doch 


— 15.272 


C2v 


3 


-22.452 


C:$ v 


-23.040 


C3 V 


4 


-30.204 


Cav 


-30.871 


Cav 


5 


-37.386 


C2v 


—38.221 


Ci 


6 


-44.891 


Cs 


-45.793 


Cs 


7 


-52.511 


C2v 


-53.456 


C2v 


8 


-60.208 


Cav 


-61.163 


Cav 


9 


-67.647 


C;ih 


-68.979 


C-sh 


10 


-75.141 


C s 


-76.398 


c 2 


11 


-82.658 


Cs 


-84.007 


Ci 


12 


-90.482 


Cav 


-91.668 


Cav 


13 


-98.672 


o h 


-99.756 


o h 


14 


-105.567 


Cs 


-106.880 


Cs 


15 


-113.132 


Ci 


-114.570 


Cs 


16 


-121.086 


C'2v 


-122.497 


C2v 


17 


-128.703 


Cav 


—129.994 


Cav 


18 


-135.761 


Ci 


-137.381 


Ci 


19 


-143.611 


Cs 


-145.255 


C2 


20 


-151.084 


Cs 


-152.667 


Cs 


21 


-158.972 


Cav 


-160.448 


Cav 


22 


-167.158 


D 4h 


-168.576 


D 4h 


23 


— 174.123 


Cs 


-175.757 


Ci 


24 


— 181.602 


Ci 


-183.331 


D 3 


25 


-189.525 


Cs 


-191.227 


Cs 


26 


-197.173 


Civ 


-198.766 


Cav 


27 


-204.809 


C'2v 


-206.602 


C2-U 


28 


-212.385 


Cs 


-214.152 


Cs 


29 


-219.969 


Cs 


-221.732 


Ci 


30 


-227.494 


Cs 


-229.315 


Cs 


31 


-235.584 


D 4h 


-237.305 




32 


-242.755 


Cs 


-244.505 


Cs 


33 


-250.295 


Ci 


-252.320 


Ci 


34 


-258.304 


Cs 


-260.187 


Cs 


35 


-265.748 


Cs 


-267.756 


Cs 
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TABLE II. The fifteen lowest energy (NaCl) 35 Cl" minima 
for the Coulomb plus Born-Meyer potential. The energies of 
the minima after reoptimization with the Welch potential 
are also given. 



Rank Energy/eV PG structure 





C+BM 


Welch 


C+BM 


Welch 












A 


1 


1 


-265.748 


-267.756 


C s 


5 


X 


5 


x 3 


B 


2 


4 


-265.741 


-267.653 


C'i 


5 


X 


5 


x 3 


C 


3 


2 


-265.739 


-267.721 


Ci 


5 


X 


5 


x 3 


D 


4 


6 


-265.685 


-267.582 


C s 


6 


X 


4 


x 3 


E 


5 


3 


-265.658 


-267.656 


Ci 


6 


X 


4 


x 3 


F 


6 


5 


-265.652 


-267.588 


Civ 


5 


X 


5 


x 3 


G 


7 


8 


-265.606 


-267.524 


c s 


6 


X 


4 


x 3 


H 


8 


7 


-265.601 


-267.573 


Ci 


5 


X 


5 


x 3 


I 


9 


12 


-265.592 


-267.476 


Cav 


8 


X 


3 


x 3 


J 


10 


25 


-265.573 


-267.417 


c s 


6 


X 


4 


x 3 


K 


11 


13 


-265.557 


-267.467 


c 2 


5 


X 


5 


x 3 


L 


12 


29 


-265.543 


-267.400 


C s 


5 


X 


4 


x 4 


M 


13 


15 


-265.542 


-267.465 


Cs 


6 


X 


4 


x 3 


N 


14 


16 


-265.519 


-267.445 


Cs 


5 


X 


5 


x 3 





15 


9 


-265.507 


-267.502 


Ci 


5 


X 


4 


x 4 



TABLE III. Details of the pathways between the lowest 
energy cuboidal nanocrystals. We give the energy difference 
between the highest energy transition state and the lowest 
energy minimum of the higher energy nanocrystal (the bar- 
rier), and the path length for the paths that minimize 
these quantities. The values of the barrier in brackets are 
for the Welch potential and were obtained by reoptimization 
of a selection of low energy pathways; an extensive search of 
the PES was not conducted with the latter potential. For 
the 5x4x4 nanocrystal we give the barriers for the minima 
L and O which are the lowest energy minima of this type 
for the Coulomb plus Born-Meyer and the Welch potentials, 
respectively. 

From To barrier/eV path length/ A 



6 


X 


4 


X 


3 




5 


X 


5 


X 


3 


0.78 


(0.69) 


33.5 


8 


X 


3 


X 


3 




5 


X 


5 


X 


3 


1.06 


(0.91) 


69.8 


5 


X 


4 


X 


4 


L 


5 


X 


5 


X 


3 


0.62 


(0.41) 


66.3 


5 


X 


4 


X 


4 


O 


5 


X 


5 


X 


3 


0.60 


(0.54) 


46.6 


8 


X 


3 


X 


3 




6 


X 


4 


X 


3 


1.06 


(0.91) 


36.3 


5 


X 


4 


X 


4 


L 


6 


X 


4 


X 


3 


0.63 


(0.44) 


57.1 


5 


X 


4 


X 


4 


O 


6 


X 


4 


X 


3 


0.56 


(0.54) 


14.4 


5 


X 


4 


X 


4 


L 


8 


X 


3 


X 


3 


1.01 


(0.83) 


92.6 


5 


X 


4 


X 


4 


O 


8 


X 


3 


X 


3 


0.97 


(0.91) 


50.7 


5 


X 


4 


X 


4 


O 


5 


X 


4 


X 


4 L 


0.60 


(0.43) 


42.8 
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